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Abstract 



Aims. Since the last phase coherent timing solution of the nearby radio-quiet isolated neutron star RX J0720.4— 3125 
six new XMM-Newton and three Chandra observations were carried out. The phase coherent timing solutions from 
previous authors were performed without restricting to a fixed energy band. However, we recently showed that the 
phase residuals are energy dependent, and thus phase coherent solutions must be computed referring always to the 
same energy band. 

Methods. We updated the phase coherent timing solution for RX J0720.4— 3125 by including the recent XMM-Newton 
EPIG-pn, MOSl, MOS2 and Chandra ACIS data in the energy range 400-1000 eV. Altogether these observations cover 
a time span of almost 10 yrs. A further timing solution was obtained including the ROSAT pointed data. In this case, 
observations cover a time span of «16 yrs. To illustrate the timing differences between the soft band (120-400 eV) and 
the hard band (400-1000 eV) a timing solution for the soft band is also presented and the results are verified using a 
Zl test. 

Results. In contrast to previous work, we obtain almost identical solutions whether or not we include the ROSAT or 
Chandra data. Thanks to the restriction to the hard band, the data points from EPIC-pn are in better agreement 
with those from MOSl, MOS2 and Chandra than in previous works. In general the phase residuals are still large and 
vary with time. In particular, the latest XMM-Newton and Chandra data show that the phase residuals have attained 
relatively large and negative values. Using this and previous timing solutions, the residuals indicate a cyclic behaviour 
with a period r; 7 — 9 yrs if the variations follow a sinusoid, or twice this value in case the residuals are modulated 
by an abs(sine) probably approaching a new minimum around MJD=55000 days (September 2009). As an alternative 
interpretation, the phase residuals can be fitted with a glitch that occured around MJD=53000 days. 

Key words, stars: individual: RX J0720.4— 3125 - stars: neutron - stars: magnetic fields - X-rays: stars 



1. Introduction 

Since the discovery of t he first isolated rad io-quiet neutron 
star RXJ1856.4-3754 (iWalter et al.||1996|) in the ROSAT 



^ All-Sky Survey only six more soft X-ray sources with sim 
ilar properties were identified. They are often referred to 
as the "Magnificent Seven" (hereafter M7; two more ob- 
jects, whi ch may be M7-lik e were found by |Rutledge et al." 



2008 and Pires et al. 2009). Unlike other isolated neutron 
stars the M7 exhibit pure thermal (black body) emission 
(with absorption fe atures in som e cases) with temperat- 



ures <100 eV (see Haberl 2007). They are conn ected to 



faint blue optical counterparts (see Kaplan 



view) with me « 25 — 28. Although these o 



2008 for a re 



in the optical bands, the optical fluxes are ss 5 — 10 times 
larger than expected by extrapolating the black body X-ray 
spectrum at low energ ies (for RBS 1774 the optical excess 
is a factor of « 30, see Zane et aL]|2008 ). 

A definite explanation of the emission properties of the 
M7 has not been put forward as yet. Their nearly Planckian 
spectrum at X-ray energies, in fact, appears difficult to re- 
concile with the predictions of standard atmospheric mod- 
els. It has been suggested that their surface layers are in 
a condensed rather than gaseous state. The phase trans- 
ition is expected to occur for low temperatures (T;S10^ K) 
and high magnetic field strengths (B^IO^"^ G; see 



TuroUa 



-: — j-l 1. et a l.|2004l and references therein and |Medin fc Lai|2006a|b 
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2007). If a thin hydrogen atmosphere covers the star, con- 
densed surface models may also explai n the optical excess 
( |Zane et al.||2007l [2004{ |Ho et al.|[2007| ). 

Among the M7, RX J0720.4-3125 (discovered by 



Haberl et al.| 1997 and identified with a faint blue optical 



star in |Motch fc Haberl||1998[ [Kulkarni fc van Kerkwijk| 
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1998| [Kaplan et al.| 2003[ [Eisenbeiss et al.| |2010[ ) has a 
unique place inasmuch it siiows significant variations in its 
surface temperature, in the equivalent width of an absorp- 
tion feature seen in t he X-ray spectrum and in the siz e of 
the emitting area (see de Vries et al.|2004 Haberl|2007 and 
Hohle et al!]|2009| hereafter H09). 

The constant spin down P of R X J0720.4-31 2 5 was 
first estimated to be « 10^^"* s/s in Zane et al. (2002) 



and an updated t iming solution using f urther observations 



was presented in Cropper et al. (|2004). The ph ase coher- 
ent timing solutions in [Kaplan fe van Kerkwijk (2005) and 



van Kerkwijk et al. (2007 hereafter KvK05 and vK07, re- 
spectively), have been computed by including either data 
from ROSAT, Chandra and XMM-Newton ("all data" solu- 
tion), or Chandra data only and all data except ROSAT 
(see vK07). Applying a constant spin-down model leads 
to large phase residuals, therefore, vK07 included a glitch 
event in one of their timing solutions which significantly 
reduced the phase residuals. 

These data sets are from various instruments with different 
energy re sponses and in differ ent data acquisition modes. 
However, Cropper et al. (2001 ) found a hardness ratio vari- 
ation and a phase shift between the flux and the hard- 
ness ratio (both folded into the 8.391 s puls e perio d) using 
XMM-Newton data. Later, |de Vries et al.| ( |2004[ ) showed 
that the energy-dependent change in the pulse profile is 
accompanied by a long-term change of the X-ray spec- 
trum an d proposed precession as a possible explanation. 
Haberl et al. ( 2004 ) confirmed the dependence of the pulse 



profiles on the energy and reported a phase lag between 
soft (120-400 eV) and hard (400-1000 eV) photons. A long 
term period of ~ 7 yrs was found in the spectral variations 
and the phase resid uals, supporting the precession model 
(Haberl et al. 20061. The phase lag exhibits a long term 
evolution, its presence appears not to depend on the tem- 
plate used to fit the X-ray pulse profile and it changed sign 
around MJD=52800 days (H09). This behaviour has to be 
taken into account and therefore requires timing solutions 
for both bands separately. Since the XMM-Newton MOS 
and the Chandra HRC and ACIS detectors are less sensit- 
ive in the soft band, we focus on a phase coherent timing 
solution for the hard band only. Since the last phase coher- 
ent timing solution in vK07, eight new XMM-Newton and 
twelve new Chandra observations were performed, thus an 
updated timing solution including the new data is required 
to distinguish between a periodic trend or a single glitch 
event. 



dius. The observations of March and September 2009 were 
analysed with SAS version 9.0.0, but following the same 
procedure as for the other observations. All observations 
are filtered using the GTI files that are generated during 
the standard data reduction process. The XMM-Newton 
observation of March 2009 was strongly affected by flares 
and the GTI files are not sufficient enough to filter out all 
time intervals with contaminated counts, i.e. for pn only 
the second half of the exposure time could be used (10.8 ks 
with « 30 X 10'^ photons in the hard band), while the MOS 
data are much less affected. This cut of the pn data is jus- 
tified posteriori, since the results (periods, phase residuals 
etc.) for pn and MOSl & M0S2 are in good agreement, as 
it will be seen later. 

The Chandra data were analysed with CIAO 4.1. For 
best corrections of the read-out times and dith er in the 
Advanced C CD Imaging Spectrometer (ACIS 
et al. [2003 ) Continuous Clocking (CC) data, we 
checked that the observations were processed with the 
Standard Data Processing (SDP) version DS 7.6.3. or later. 
The coordinate accuracy for all ACIS-CC observations is 
better than 0.5", i.e. one pixel. We applied the Charge 
Transfer Inefficiency (CTI) correction and use the source 
photons from a rectangular region of 1" x 1" size covering 
the brightest pixels (located on chip seven for all AC IS-CC 
observation s) . For High-Resolution Camera (HRC, iKraft 



Garmire 
first 



et al.||1997p data we used photons from zeroth order 
cular 



cir- 

within 



region of 2" radius) and both first orders 
the standard LETG spectral extraction windows (but lim- 
ited to 10 — 60 A and the brightest 1" pixel stripes only). 
Finally we applied the barycentric correction using ELxbary 
for the HRC and ACIS-CC data. 

The ROSAT data were taken froml Cropper et al.| ( [2004| 
whom we refer to for details of the observations and the 
data reduction. 

Assuming that no significant variations in the period 
of RX J0720.4— 3125 occur within a few days, we merged 
several close Chandra observations with the same instru- 
mental setup to reduce the scatter in the phase folded light 
curves. After this merging we obtain 70 data sets: 16 from 
XMM-Newton EPIC-pn and 36 from EPIC-M0S1/M0S2 
(the XMM-Newton observations rev. 0622 and 0711 from 
May and October 2003, respectively, have four MOS data 
sets each), 6 ROSAT data sets (pointings only) and 12 (out 
of 29) data sets from Chandra observations. All observa- 
tions arc listed in Tab. [l] first column. 



2. Observations and Data reduction 



3. The selection of the hard band (400-1000e\/) 



We use EPIC (pn [Striider et aL] [200r, MOSl and M0S2, 
Turner et al.|2001[) data from all 14 XMM-Newton observa- 
tions performed between May 2000 and November 2007 (for 
a summary of the instrumental setups see H09) and include 
the most recent XMM-Newton observations of March 2009 
(rev. 1700, performed in full frame mode with the thin fil- 
ter for pn and both MOS) and September 2009 (rev. 1792, 
performed in full frame mode with the thin filter for pn and 
in small window mode with thin filter for both MOS). We 
analysed the XMM-Newton data with the standard XMM- 
Newton Science Analysis System (SAS) version 7.1.0. ap- 
plying barycentric correction. We use single and double 
events for both MOS and the pn detector, that are collec- 
ted within a circular region of 30" (pn) and 8" (MOS) ra- 



in previous work ( [Haberl et al.j2004[ [20061 ^'^'^ H09) a vari- 
able phase lag between soft (120-400 eV) and hard (400- 
1000 eV) photons of RX J0720.4-3125 was reported and 
discussed. Depending on the different detectors and acquis- 
ition modes, the data contains different fractions of photons 
from the two bands that lead to a systematic scatter in the 
phase residuals. This phenomenon has to be taken into ac- 
count and requires a different treatment as for the other 
M7 if an adequate timing solution has to be found. 

In this section we first discuss the cross calibration issues 
for the different instruments used for the observations of 
RX J0720.4— 3125, then derive an empirical energy band 
that minimises the scatter in the residuals and discuss the 
results. 
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Figure 1: The accordance of the phase residuals of the 
hard band (cut off energy- 1000 eV) for XMM-Newton 
EPIC-pn, EPIC-MOSl & M0S2 and Chandra HRC of 
RX J0720.4— 3125 for different values of the energy band 
cut. Data sets taken with different observatories were se- 
lected in such a way that they are as close in time as pos- 
sible. The upper panel shows the sum of the absolute val- 
ues of the phase differences (phase gaps) derived from the 
different instruments with respect to each other, while the 
middle panel shows the median of these phase gaps and the 
corresponding x^/d.o.f. is presented in the lower panel. 
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Figure 2: The same as in Fig [T| but using the soft band 
(120 eV— cut off energy) for the analysis. 
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In |Haberl et al.| (|2004) |2006| the photons were divided 
into a soft (120-400 eV)andahard band (400-1000 eV) for 
the EPIC-pn detector to illustrate the phase shift between 
these two bands. The effective area of the MOS detectors 
deviates much from that of the pn detector (particularly 
in the soft band) and the spectral res olution of the MP S 
detectors is ~ 100 eV in the soft band ( jTurner et al.l200lL 
This resolution is much worse than for gratings of A.MM- 
Newton or Chandra. In addition, the CCD detectors suffer 
under the contamination of hard photons in the soft band 
due to redistribution. The MOS and pn observations were 
executed in different filters: thin, medium and thick. The 
thick filter significantly suppresses the penetration of soft 
photons compared to e.g. the thin filter, that also influences 
the timing properties. 

The phase residuals from the XMM-Newton observa- 
tions should be comparable to those of the Chandra ob- 
servations, taking different energy response, spectral resol- 
ution and the phase shift between hard and soft photons 
into account (softer photons are less prominent in the MOS, 
HRC and ACIS data than in the pn, influencin g the resid- 
uals as shown in H09, Fig. 1 and |Haberl|2007| Fig. 7). The 
HRC observations have a low number oi photons and the 
energy information of the photons in the zeroth order is 
lost. Therefore we included all photons in the HRC data. If 
we would exclude the zeroth order photons, ^ 40% of the 
counts were lost and if only the hard band (first orders) 
would be used, ~ 80% of the counts were lost, resulting 
too low statistics. For similar reasons (low statistics, low 
energy resolution) we also used ROSAT data from a wider 
energy band: channels 10 to 90 (100-900 eV) for the Position 
Sensitive Proportional Counter (PSPC) and channels one 



to eig ht for the High Resolution Imager (HRI, David et al 
1996f . 

The energy band minimising the residuals can be 
found using those Chandra observations performed close 
to XMM-Newton observations: The HRC observation 5582 
(June 1, 2005) is close to the XMM-Newton observation rev. 
0986 (April 28, 2005 ), the HRC observations 6369 & 7177 
(Oct 8/9, 2005) are close to the XMM-Newton observation 
rev. 1060 (Sep 23, 2005), the HRC observations 7243-7245 
(Dec 14-17, 2005) are close to the XMM-Newton observa- 
tion rev. 1086 (Nov 11/12, 2005), the HRC observations 
10861 & 10700 (Jan 20 & Feb 14, 2009) are close to the 
XMM-Newton observation rev. 1700 (Mar 21, 2009) and 
the HRC observation 10701 (Sep 11, 2009) is close to the 
XMM-Newton observation rev. 1792 (Sep 22, 2009) see also 
Tab.jlj column 6 and 7. Unfortunately no ACIS observation 
is close in time to these observations. 

Assuming that the timing properties of 
RX J0720.4— 3125 do not change significantly within 
a time span of a month, we compare the timing residuals 
(using the "all data" timing solution in vK07) of these 
observations directly changing the band cut from 200 eV 
to 800 eV in steps of 10 eV. For example, if the band 
cut is at 600 eV, we have two bands of 120-600 eV and 
600-1000 eV. Due to the limited energy resolution, these 
band cuts could not be applied for the HRC observations. 

We then calculate the phase residuals by fitting a sine 
for the phase folded light curves (see vK07) in 12 phase bins 
for pn, MOSl and M0S2 and HRC for each band cut and 
calculate the sum of the absolute values and the median of 
the phase gaps of the phase residuals of the harder band 
from the different detectors. Having the phase residuals and 



their errors, we also calculate x^/d.o.f showing the degree 
of accordance of the results from the different detectors. In 
all three cases, we derive a minimum (i.e. best agreement) 
if we use the band cut between 300 eV and 400 eV (for 
consistency with previous work, |Haberl et al.||2004l |2006[ 
and H09, we set the cut between the two bands at 400 eV), 
i.e. having two bands of 120-400 eV and 400-1000 eV and 
using the harder band (see Figure IT]). For band cuts of lower 
energies the different detector responses for the soft photons 
cause larger phase gaps, because of the phase shift between 
hard and soft photons and the rough energy resolution of 
the EPIC detectors, while for band cuts at higher energies 
the number of photons decreases and the statistics worsens. 

If the same is done, but keeping the soft band photons 
(see Fig. [2]) , the best agreement of the different instruments 
is achieved if no band cut is applied. But even the lowest 
X^/d.o.f value is more than two times larger as in the case 
when photons between 400-1000 eV are used. 

Apart from the intrinsic properties of RX J0720.4— 3125 
(the spectrum did not change much since rev. 0986, see 
H09), the choice of this band cut for the hard photons at 
least reduces instrumental discrepancies of the different de- 
tectors. 



4. Timing procedure 

For comparison wit h the results of KvK 05 (obtained us- 



ing the Z^ method, Buccheri et al. 1983), we start taking 



as reference period that from the ACIS-CC observations 
4666-4669 combined with the HRC observation 5305, cov- 
ering a total time span of 52 days. We determine the period 
from the peak in the Z^ periodogram (from now on always 
n = 1) and estimated the la errors f ollowing Kaplan et aLj 
( |2002 ) using the equation derived in Ransom et al. (2002| T^ 
The period we have taken corresponds to the center of the 
la region in the periodogram (since the peak is symmetric 
close to its maximum). 

Using this dataset produces aliases in the periodogram 
caused by the time gaps between the observations. We ob- 
tain three peaks with Zj w 140: Pi = 8.39036664(53) s, 
P2 = 8.39111600(50) s and P3 = 8.39188114(53) s. P2 yields 
formally the largest Zf value and is consistent with the ref- 
erence period P = 8.39111590(50) s from KvK05 (numbers 
in parentheses indicate la errors using the error estimation 
in [Ransom et al. 2002 ) . During the timing procedure we 
found F2 being the correct period since starting with Pi or 
P3 as initial period, the procedure does not converge. 

As a first step we determined the period of each obser- 
vation separately using the Zf method. The period of the 
Chandra HRC observation 7251 shows a large discrepancy 
with respect to the others, probably caused by the small 
number of photons giving not enough statistics for a reli- 
able period determination. Therefore we excluded this ob- 
servation from our analysis. The periods obtained from the 
ROSAT data with a low count number differ sensibly from 
the reference period and deviate significantly from those 



periods listed in Cropper et alT] ( [200 4) and are highly un- 
certain, thus are also excluded for the first phase coherent 
timing solution. No periodic signal was found in the ROSAT 
all sky survey, i.e. this observation is not used at all and 
not listed in Tab. [T] 

Most of the individual periods with more than 10^ 
photons are cons istent with the sing le spin-down model 
within la errors (Ransom et al. 2002) and all of them are 
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Figure 3: Phase binned (40 bins) light curve (400-1000 eV) 
of the XMM-Newton EPIC-pn observation revolution 1086 
(November 2005) folded during the timing procedure (de- 
scribed in the text). The solid line represents the combin- 
ation of three Fourier harmonics chosen for the determina- 
tion of the phase. The phase shift differs slightly from that 
from the sine fit (dotted line). The dashed line shows, that 
a fit with the second harmonic (cosine) is not always suffi- 
cient. The error bars denote Poissonian errors. 



consistent within 2ct errors, both for the hard and the soft 
band. Note that the periods for Chandra HRC and ROSAT 
data are the same for both bands since we applied no en- 
ergy selection for these observations, as mentioned before. 
All individual periods are listed in Tab. [2] column 4 and 5. 

Our timing procedure follows that in KvK05 with the 
difference that we calculate the times of arrival (TOA, the 
time of maximum light closest to the middle of the obser- 
vation, see KvK05) from the included observations during 
the iterations of the timing procedure. The errors of / and 
P are evaluated after each step of the iteration: in the case 
of / the error is calculated from the different / values satis- 
fying the condition to minimise the phase residuals within 
their la errors. This restricts the possible values for / if 
more observations are added, i.e. the error on / decreases 
with growing time T from the reference point. The error 
of P, eP, automatically decreases with growing time span 
T following eP ~ 1/T (e.g. see the general approach in 
Kovacs|p8l] ). 
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Figure 4: The Fourier coefficients (cosl=0 by definition) of 
the fit to the phase folded light curves (40 bins) of the 
EPIC-pn observations. The values of the coefficients are di- 
vided by the constant offset Aq, since the observations have 
different numbers of photons. The first panel shows the cor- 
responding values of x^/d-o.f.. All observations have less 
than 1% background counts and are filtered using standard 
GTI files. In case of rev. 1700 {^ 43% high background) an 
additional filtering is applied (see Sec. 2) leading to the de- 
viating point near MJD=54900 days. All error bars denote 
Ict confidence level. 



the Fourier serie^as follows: The Fourier series for the light 
curves is given by 



F(0) = ^ AkSin(k0) + Bkcos(k0) 



(1) 



k=0 



In KvK05, vK07 and II09 the phase binned light curves 
were fitted with a sinusoid to derive the phase of the max- 
imum, although the pulse p rofiles are not alway s best rep- 
resented by a sine function ( Haberl et al.|[2004 |. We fitted 
the light curves using a combination of three harmonics of of fit are scaled to \Jx^ jdoj . 



where k = corresponds to the constant offset. We 
choose to fit the light curves with the sum of three terms, 
not to introduce too many degrees of freedom. To sim- 
plify the notation, in the following we denote each term 
just as sink or cosk, e.g. sinl stands for Aisin(10), cosl for 
Bicos(10). Since the contribution of cosl turned out to be 



^ The light curves are fitted using the MATLAB® (version 
R2008a) internal function fit with default adjustments. Errors 
derived from this function and the errors derived from the results 
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Figur e 5: The principal components (see Zane & Turolla 
2006 for details) of the 16 EPIC-pn observations of 
TDU0720.4-3125. 



negligible (since we forced the phase to be zero at the refer- 
ence time, and a cosine term would imply a phase shift), our 
sum starts always with sinl, while the remaining two terms 
can be sin2-|-sin3, cos2-|-sin3 or cos2+sin2, and so on. To 
keep the number of free parameters small, we do not use 
terms with k > 3 and use only three terms in total. The 
combination that fits the individual light curve with the 
lowest value of x^/d.o.f. was used for the phase determina- 
tion. This leads to an improvement of the x^/d.o.f. for the 
light curves (see example in Fig. p| and to a better determ- 
ination of the phase shift. Our approach is justified since 
in most of the cases some Fourier coefficients are negligible 
(see the coefficients from the phase folded light curves of 
the EPIC-pn observations in Fig. Has an example). 

The analysis of the Fourier coefficients shows that the 
pulse profile chan ges in time, as a l ready reported in lHaberl] 



et al. (2006) and de Vries et al. (2004 their Fig. 2), and 



the evolution continued also after the last observation con 



sidered by [Haberl et ah] ( |2006[ MJD=53700 days). In order 
to check that the light curve evolution is independent of the 



adopted pu lse profile template (either a sine as in Haberl 
|et al.|[2i306| or a truncated Fourier expansion, as in the 
present work), we performed a principal component ana- 
lysis (PCA, see e.g. Zane fc Turolla|[2006 for more details) 
of the same 16 EPIC-pn light curves we used before, again 



binned in 40 phase intervals. Starting from the original vari- 
ables (the 40 values of the counts at the different phases), 
the PCA computes a new set of variables (the principal 
components, pes) which are a linear combination of the old 
ones and are ordered in such a way that the first pc accounts 
for the largest variance of the data, the second the second 
largest, and so on. We find that the first four pes are actu- 
ally responsible for ^ 92% of the variance, pulse profiles can 
be effectively classified in terms of only a limited number 
of pes which embody their main characteristics. The time 
evolution of the first four pes is shown in Fig. [5 The signi- 
ficant changes in both sets of coefficients (pc and Fourier) is 
a strong evidence of a genuine variation of the pulse profile 
in the hard band in time. The second last XMM-Newton 
observation (rev. 1700) that was contaminated by strong 
flares de viates from the trends i n Fig. |4] and Fig. [5| As dis- 



cussed in Zane & Turolla (20061, the first pc (Zi) is related 
to the pulse amplitude, Z2 provide a measure of the phase 
interval of the light maximum and Z3 reflects the li ght curve 
parity with respect to the half period (see Fig. 4 of Zane & 
[Tur olla 2006, for the dependence on the phase of the first 
four coefficients Vik, with i — 1,...,4, used to calculate 
the PCs from the original variables). Although a detailed 
analysis will not be attempted here. Fig. [5] shows that the 
largest changes occur in Zi and Z2 and are then associated 
to the amplitude and the position of the maximum of the 
pulses. 



5. Results 

We performed the phase coherent timing by binning the 
light curves into 10, 12, 16, 18 and 20 phase bins. A large 
number of bins leads to large scatter in the Chandra light 
curves due to the small count rate compared to XMM- 
Newton while too few bins result in an insufficient time 
resolutiorjj We derived the best timing solution for all ob- 
servations (excluding ROSAT data and the HRC observa- 
tion 7251) for 12 phase bins in order to minimise the phase 
residuals. 

If we fit the phase folded light curves with a pure sinus- 
oid, we obtain P = 8.3911153307(22) s, / = -9.933(52) x 
10^"'^^ Hz/s for the phase coherent timing solution. This 
corresponds to y^ /d.o.f. — 12.7 for the timing solution 
and for the light curve fits we obtain y^ /d.o.f. = 2.29 on 
average. Fitting the light curves with a variable combin- 
ation of the three Fourier harmonics, the result does not 
change significantly. We obtain P — 8.3911153362(39) s 
and / = -9.946(74) x lO^^*^ Hz/s, but the phase resid- 
uals have smaller errors and therefore the timing solution 
corresponds to x^ jd.o.f. — 47. For the light curve fits we 
obtain ^ jd.o.j. — 1.23 on average. 

Although the periods d etermined from the R OSAT data 
(see Zane et al. 2002 and Cropper et al.| 2004 for details) 
differ significantly to the reference period P2, we include 
these data sets, except the All-Sky Survey (where we find 
no period), for a further timing solution. The inclusion 
of the ROSAT data extends the time span of observa- 
tions from ^9.6 yrs to «16.5 yrs and enlarges the data 



■^ Note that the time resolution for the MOS detectors is 0.3 s 
(small window), 0.9 s (large window) and 2.6 s (full frame), 
i.e. less than one phase bin. This is compensated by the large 
number of photons. KvK05 used 16 phase bins (including the 
MOS detectors). 
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I values in the P — f plane derived from the XMM-Newton, ROSAT and Chandra observations (except 



Figure 6: The Z? 

the HRC observation 7251) of RX J0720.4— 3125 compared to the results of the phase coherent timing solutions in this 

work. The peak is located at P = 8.391115309(14) s and / = -9.992(61) x IQ-^^ Hz/s. The three different timing 

solutions derived without the Chandra data (open circle), without ROSAT data (star) and using XMM-Newton, ROSAT 

and Chandra observations (dot) are not resolved on this scale (red box), but become visible enlarging the red box (shown 

in the upper left panel); the values are listed in Tab 

solutions are derived from the hard band (400-1000 eV). 



The errors of the Z\ solution are scaled to \/x^~[dof . All timing 
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set from 64 to 70 observations. Again, we fitted the light 
curves with a variable combination of the three Fourier 
harmonics and obtain P = 8.3911153336(22) s and / = 
-9.961(67) X 10"^^ Hz/s for 12 phase bins. This corres- 
ponds to x^ /d-o.f. — 45 for the timing solution and for the 
light curve fits we obtain x^/d-o.f. = 1.17 on average. 

For a third timing solution we excluded the Chandra 
data, but using the ROSAT data (58 data sets out of 70) 
and the same reference time and period as for the other 
two timing solutions. Fitting the light curves with a vari- 
able combination of the three Fourier harmonics, we obtain 
P = 8.3911153310(22) s and / = -9.940(71) x lO^^^ Hz/s 
for 12 phase bins that corresponds to y^ /d.o.f. — 47 for 
the timing solution and for the light curve fits we obtain 
X^ /d.o.f. = 1.21 on average. 

To verify the correctness of our timing solutions we 
apply the Zj method to the same combination of obser- 
vations as for the phase coherent timing solutions (see 
Fig. [6]). Using all 70 data sets, we found a maximum at 
P = 8.391115309(14) s and / = -9.992(61) x lO^^^ Hz/s 
with Zf = 8263. The errors for P and / are obtained from 
the equation in Ransom et al. (2002), whereas the values 
correspond to the maximum peak in the center of the la 
confidence region. In this region the peak is symmetric. 
Finally, the errors are scaled to \/x^/dof. 

The results of the three phase coherent timing solutions 
and from the Z^ method are summarised in Tab. pi Note 
that our results are only slightly different from previous 
solutions (but scatter less), i.e. the restriction of our in- 
vestigations to the hard band reduces the phase residuals 
and influences the shape of their evolution in time, but does 
not affect basic parameters like the spin-down age or mag- 
netic field strength of RX J0720.4-3125 significantly. For 
completeness the results from the coherent " all data" solu- 
tion and the Zf method for the soft photons are listed in 
Tab. |3] too, illustrating the difference of the timing solu- 
tions in the different bands. The corresponding TOAs from 
the final solutions are listed in Tab. [2] last two columns. 
The errors of the TOAs are derived from the formal fit er- 
rors of the phase binned light curves by fitting them to the 
combination of Fourier harmonics as explained in Section 
4. 



6. Discussion 

We present new updated phase coherent timing solutions 
for the isolated radio-quiet X-ray pulsar RX J0720.4— 3125 
including the most recent XMM-Newton and Chandra ob- 
servations, restricting our analysis to the hard energy band 
(400-1000 eV), except for Chandra HRC and ROSAT data. 
The new solutions were obtained with and without the in- 
clusion of the ROSAT and Chandra data and differ slightly 
from previous ones (see KvK05, and vK07 and Tab. l3|. 
We improved the phase determination and the modeling of 
the phase folded light curves by fitting a variable combin- 
ation of the three Fourier harmonics with the sine as the 
leading term, instead of a pure sinusoid. We checked the 
phase coherent timing solutions applying the Zj method 
to the observations in the P/f plane. All timing solutions 
correspond to a -y^ /d.o.f. « 50, i.e. the phase residuals 
are still large and a timing solution based only on / and 
P provides an inadequate model for the spin behaviour of 
RXJ0720.4-3125. 



-0.1E- 
4.9 




sine: xl, = 5.29 P,^^ = (8.91 ± 0.21 ) yrs 
abs(sine): xj,,, = 6.52 P ^^ = (17.48 ± 0.24) yrs 



5.1 5.2 5.3 5.4 5.5 5.6 

MJD[days] ^.,0^ 



Figure 7: Phase residuals after applying the timing solutions 
derived without the ROSAT observations and without the 
Chandra HRC observation 7251 (marked as a filled circle) 
from RX J0720.4— 3125 applied to all observations listed in 
Tab. [l] The best fits with a sine and an abs(sine) model 
(only the observations used in the current timing solution 
are used for the fits) are also shown. Squares mark resid- 
uals derived from XMM-Newton EPIC-MOSl and M0S2 
observations, dots EPIC-pn observations, circles Chandra 
HRC and ACIS-CC data and stars ROSAT data. Ah error 
bars denote Icr. 



Due to the restriction to the hard band for most of 
the data, the phase residuals from EPIC-pn do not devi- 
ate systematically to those from EPIC-MOSl and M0S2 
and Chandra unlike those derived without energy restric- 
tion (see Fig. 1 in H09 for comparison and Fig. IT] in this 
work) and most observations are in good agreement to each 
other. However, the residual of the Chandra observation 
7251 is still not consistent (lowest number of counts and bad 
statistics) to the other phase residuals and the phase resid- 
uals from the XMM-Newton EPIC-pn and EPIC-MOSl & 
M0S2 observations 0622 scatter much. The observations of 
rev 0622 were all performed in different filters: thick, me- 
dium and thin, having different transparency for the softer 
photons (causing deviations in the phase as discussed be- 
fore) in the hard band. 

Including the most recent Chandra and XMM-Newton 
observations, the ROSAT data become less important for 
the timing solution. This explains why our three phase co- 
herent timing solutions are almost identical and scatter less 
than the phase coherent timing solutions in KvK05 and 
vK07 (without a glitch). 

Applying the timing solutions to the observations of 
RX J0720.4— 3125 we still obtain large residuals. Due to 
the small number of counts, the ROSAT data are ambigu- 
ous in determining the correct phase, the periods differ sig- 
nificantly from the periods of the other observations and 
the ROSAT periodograms are noisy. However, the inclu- 
sion of the ROSAT data extends the observed time span 
from 9.6 to 16 yrs. Therefore, we present the phase re- 
siduals from all observations, but derived from the tim- 
ing solution without ROSAT data (hard band) in Fig. [t] 
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Table 3: P and / derived from different methods in this paper, compared to previous results. 



solution 


period-8.39flf5 


./ 


rms 


d.o.f 


xVd.o.f 




xfO"'' s] 


xlO"" [Hz/s] 


[s] 








this work, 


hard band 








"all data" 


3.336(22) 


-9.961(67) 


0.62 


70-3 


45 


after fitting a sine 






0.29 


70-5 


7.8 


after fitting an abs(sine) 






0.31 


70-5 


8.6 


without ROSAT 


3.362(39) 


-9.946(74) 


0.60 


64-3 


47 


after fitting a sine 






0.19 


64-5 


5.3 


after fitting an abs(sine) 






0.21 


64-5 


6.5 


without Chandra 


3.310(22) 


-9.940(71) 


0.61 


58-3 


47 


after fitting a sine 






0.31 


58-5 


9.7 


after fitting an abs(sine) 






0.33 


58-5 


8.7 


Z^ "all data" 


3.09(14) 


-9.992(61) 


0.73 


70-3 


64 


Zf without ROSAT 


2.96(13) 


-10.047(34) 


0.58 


64-3 


46 


Zi without Chandra 


3.09(13) 


-9.980(36) 


0.72 


58-3 


51 




this work. 


soft band 








"all data" 


3.429(22) 


-9.956(72) 


0.50 


70-3 


40 


Zl "all data" 


3.31(62) 


-9.959(17) 


0.50 


70-3 


39 


previous work (appli 


3d to the hard band) 






vK07 ("all data") 


2.670(84) 


-9.88(13) 


0.97 


70-3 


81 


vK07 (without ROSAT) 


2.846(77) 


-9.74(04) 


1.30 


64-3 


207 


KvK05 ("all data") 


3.20(13) 


-9.918(15) 


0.64 


70-3 


56 


KvK05 (Chandra) 


3.05(16) 


-9.97(06) 


0.64 


12-3 


57 



Note: For the timing solution in this work we always excluded the Chandra HRC observation 7251. Since the phase residuals seem 
to follow a periodic patt ern we fitted a sine and an abs(sine), see Sec. 6. All errors correspond to la confidence (for the Zf solution 
see Ransom et al. 12002 



with errors scaled to Wx'^/d-of, for the phase coherent solutions, see Sec. 5). 



shown with an error weighted sine and an abs(sine) fit (see 
also H09). The phase residuals reached a minimum around 
MJD=52800 days and continue towards negative values, 
maybe a new minimum around MJD=55000 days; but this 
has to be confirmed by new observations finding a possible 
turning point. There is no clear evidence for a minimum in 
the past (probably MJD=50000 days including the ROSAT 
data), but the ROSAT data are less reliable. 

Fitting a sine or an abs(sine) to the phase residuals re- 
duces the x^/d-o.f. from w 50 to w 6 — 7 with a period 
Pros ~ 8.9 yraj (i.e. the observations cover two cycles) or 
twice this period fitting an abs(sine), respectively. However, 
the -y^ /d.o.f . is formally still unacceptable. A sinusoidal be- 
haviour of the phase re siduals could be expla ined by pre- 
cession, as s uggested in de Vries et al. (2004) and Haberl 
et al. (2006). Precession would cause an advanced (the NS 



precesses towards the observer) and a retarded (the NS 
precesses backwards with respect to the ob server) signal, 
where as the residuals would follow a sine (Nelson et al. 
19901. During precession the observer would see different 



parts of the surface that would lead to changing spectral 
properties, such as changes in temperature and size of the 
emitting surface. If two, roughly antipodal, hot spots would 
both contribute to the X-ray emission, the phase residuals 
caused from both hot spots would be shifted by Pres/2: 



^ In H09 long term periods of Prcs ~ 5.5 — 7.5 yrs were found 
for the sine fit. This is significantly less than found in this work. 
However, the timing solutions used in H09 were not derived from 
the inclusion of the new data but were adopted from KvK05 and 
vK07, also required the inclusion of an additional slope (not 
necessary in this work). 



if one hot spot precesses towards the observer and ap- 
pears, the second hot spot moves backwards with respect 
to the observer and disappears. Since the timing is not 
sensitive to which of the hot spots is causing the residuals, 
the phase residuals qualitatively may have the shape of an 
abs(sine), having twice the period of the corresponding sine 
(two peaks from two different hot spots). The presence of 
the larger period is supported by th e variable pha s e lags 



between hard and soft photons, see Haberl et al. (2004 
2006 and H09). However, it has to be shown whether this 



scenario can explain the spectral behaviour in H09 too. 
Formally, then also the spectral changes undergo a periodic 
behaviour having P^es ~ 2 x 8.9 yrs. A long term period of 
Pres ~ 17 — 18 yrs could explain why the spectral changes 
do not show a periodic trend (H09) yet, since the time span 
of the XMM-Newton observations would cover less than one 
cycle. The new turning point of the phase residuals would 
be at MJD = 553871^^^ days (first half of July 2010) or at 
MJD = 552401^ days (first half of February 2010) for the 
sine or the abs(sine) fit, respectively (see Fig. [T]). 

Alternatively, vK07 proposed a glitch around 
MJD=53000 days explaining the spectral and temporal 
changes. As shown in H09, the temperature increased 
at this time, followed by a slow decrease, as expected 
during and after a glitch event. Therefore vK07 presented 
a "glitch solution" to explain also the timing behaviour of 
RX J0720.4— 3125. This glitch solution was not applied cor- 
rectly in H09. Indeed, the glitch solution presented in vK07 
minimises the phase residuals for the observations avail- 
able at that time (MJD=53500 days), i.e. the statement 
(in H09) that it produces large residuals is not correct. 
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However, adding new observations after MJD=:53500 days 
(not available in vK07) and applying the glitch solution 
correctly, the phase residuals grow to larger values. This 
can be shown even with the observations available in H09 
(MJD=54421 days) and is even more pronounced using the 
recent data (until MJD=55100 days), i.e. the conclusions 
in H09 remains the same. 

To reduce the residuals, we introduce a new additional 
parameter with the physical meaning of a post-glitch cor- 
rection for /. This can be done by fitting a parabola with 
respect to the time distance to the proposed glitch time, 
tg = 52866 days, in vK07 ("aU data" solution). The er- 
ror weighted fit for the hard band including all data gives 
fc = —1.11(10) X 10^^^ Hz/s and significantly reduces the 
phase residuals {x^/d.o.f. = 2.8, rms=0.31 s). vK07 ob- 
tained fvK07 = —1.04(3) X lO"^'^ Hz/s after tg, the new 
value would be fnew — .fvK07 + fc, i-e. a modification of the 
glitch solution in vK07. Due to this correction the glitch 
hypothesis is still a competing model. 

The phase residuals of RX J0720.4— 3125 after applying 
the glitch solution in vK07 and the post-glitch correction 
(both for the hard band) are shown in Fig. M 

Recentl y, an analysis of tim ing irregularities of 366 ra- 
dio pulsars ( Hobbs et al.|2010 ) showed that quasi— periodic 
structures in the timing residuals are quiet usual - at 
least for radio pulsars. These quasi— periodic structures 
are dominant for young pulsars (age< 10^ yrs, note that 
RX J0720.4— 3125 is probably much younger than implied 
by its char acteristic age, see |Tetzlaff et al.|[20Tm |Kaplan| 



et al. 



2007 1 . In some cases the amplitude of the timing re- 



siduals are comparable to those of RX J0720.4— 3125 and 
the timing residuals show "periods" of a few years. These 
radio pulsars were studied over decades and the data points 
are much denser in time than in the case of X-ray pulsars, 
thus the timing residuals of RX J0720.4— 3125 may also fol- 
low such trends. In such cases (as the authors argue in 
Hobbs et al.|2010 1 the timing residuals caused due to " slow 



glitches" may not be a different phenomenon to that caus- 
ing the timing irregularities. 

The behaviour of RX J0720.4-3125 is stiU not under- 
stood. The monitoring of RX J0720.4-3125 using XMM- 
Newton and Chandra is still ongoing and will help to bring 
us closer to understand the reason for its behaviour. 
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Table 1. All observations of RX J0720.4— 3125 listed in chronological order. 



Obs. Id. 


Instrument/setup 


Counts 




Counts 


MJD 


Start Date 


Effective 






(120-400 cV) (400-1000 cV) 






exposure 












[days 




[ks 


rp300338n00 


PSPC 




5374 




49258 


1993 Sep 27 


3.22 


rh300508n00 


HRI 




1259 




50199 


1996 Apr 25 


3.13 


rhlSOlOOnOO 


HRI 




1197 




50211 


1996 May 7 


3.57 


rh300508n01 


HRI 




493 




50354 


1996 Sep 27 


1.41 


rh400884n00 


HRI 




13381 




50391 


1996 Nov 3 


33.57 


h400944n00 


HRI 




3054 




50924 


1998 Apr 20 


3.57 


368^ 


HRC-S/LETG 




2722 




51575 


2000 Feb 1 


5.40 


7451 


HRC-S/LETG 




9392 




51576 


2000 Feb 2 


26.26 


369^ 


HRC-S/LETG 




2660 




51579 


2000 Feb 4 


6.12 


0078 S3 


EPIC-pn/FF thin 


241783 




115287 


51677 


2000 May 13 


58.60 


SI 


EPIC-MOSl/FF thin 


69623 




43977 


51677 


2000 May 13 


61.98 


S2 


EPIC-M0S2/SW thin 


68010 




40818 


51677 


2000 May 13 


61.99 


0175 S3 


EPIC-pn/FF mod 


103491 




55743 


51870 


2000 Nov 21 


25.65 


S7 


EPIC-MOSl/LW mcd 


15765 




11662 


51870 


2000 Nov 21 


18.00 


27742 


ACIS-CC 


13455 




15648 


52248 


2001 Dec 4 


15.01 


27732 


ACIS-CC 


9452 




11646 


52248 


2001 Dec 5 


10.61 


27712 


ACIS-CC 


876 




1553 


52250 


2001 Dec 6 


1.86 


27722 


ACIS-CC 


3537 




4487 


52250 


2001 Dec 6 


4.05 


0533 S3 


EPIC-pn/FF thin 


137951 




72150 


52585 


2002 Nov 6 


28.38 


SI 


EPIC-MOSl/FF thin 


27987 




21892 


52585 


2002 Nov 6 


29.99 


S2 


EPIC-M0S2/FF thin 


29169 




22192 


52585 


2002 Nov 6 


29.99 


0534 S3 


EPIC-pn/FF thin 


145057 




76095 


52587 


2002 Nov 8 


30.18 


SI 


EPIC-MOSl/FF thin 


28168 




22112 


52587 


2002 Nov 8 


31.80 


S2 


EPIC-M0S2/FF thin 


28954 




22768 


52587 


2002 Nov 8 


31.79 


0622 U2 


EPIC-pn/SW thick 


123199 




94655 


52762 


2003 May 2 


72.79 


S5 


EPIC-MOSl/FF med 


21204 




18741 


52762 


2003 May 2 


29.98 


S3 


EPIC-MOSl/FF thin 


28433 




23754 


52762 


2003 May 2 


33.59 


S6 


EPIC-M0S2/FF med 


21988 




19650 


52762 


2003 May 2 


29.99 


S4 


EPIC-M0S2/FF thin 


30901 




24413 


52762 


2003 May 2 


33.59 


0711 S8 


EPIC-pn/SW mcd 


83151 




60557 


52940 


2003 Oct 27 


24.90 


S3 


EPIC-MOSl/FF thin 


10841 




11937 


52940 


2003 Oct 27 


13.90 


S5 


EPIC-MOSl/LW thin 


10547 




11769 


52940 


2003 Oct 27 


15.71 


S4 


EPIC-M0S2/FF thin 


11649 




12392 


52940 


2003 Oct 27 


13.91 


S6 


EPIC-M0S2/LW thin 


11943 




12382 


52940 


2003 Oct 27 


15.71 


46663 


ACIS-CC 


5606 




13094 


53010 


2004 Jan 6 


10.12 


46673 


ACIS-CC 


2678 




5835 


53011 


2004 Jan 7 


4.79 


4668^ 


ACIS-CC 


2121 




4918 


53016 


2004 Jan 11 


5.16 


4669^ 


ACIS-CC 


1890 




4277 


53023 


2004 Jan 19 


5.22 


5305^ 


HRC-S/LETG 




13275 




53062 


2004 Feb 27 


35.70 


0815 SI 


EPIC-pn/FF thin 


130683 




93280 


53147 


2004 May 22 


41.30 


S2 


EPIC-MOSl/FF thin 


35853 




41178 


53147 


2004 May 22 


45.21 


S3 


EPIC-M0S2/FF thin 


35915 




42318 


53147 


2004 May 22 


45.22 


4670^ 


ACIS-CC 


5723 




12329 


53221 


2004 Aug 3 


10.13 


4671^ 


ACIS-CC 


2778 




6318 


53223 


2004 Aug 5 


5.15 


4672^ 


ACIS-CC 


2483 




5789 


53227 


2004 Aug 9 


5.12 


4673^ 


ACIS-CC 


2527 




6034 


53244 


2004 Aug 23 


5.13 


5581^ 


HRC-S/LETG 




30998 




53393 


2005 Jan 23 


68.20 


0986 S3 


EPIC-pn/FF thin 


181754 




131146 


53488 


2005 Apr 28 


51.43 


SI 


EPIC-MOSl/SW thin 


32400 




36475 


53488 


2005 Apr 28 


53.05 


S2 


EPIC-M0S2/SW thin 


34456 




38266 


53488 


2005 Apr 28 


53.06 


5582^ 


HRC-S/LETG 




35777 




53523 


2005 Jun 1 


70.17 


63648 


HRC-S/LETG 




22210 




53610 


2005 Aug 27 


38.87 


1060 S3 


EPIC-pn/FF thin 


175481 




120913 


53636 


2005 Sep 23 


51.14 


SI 


EPIC-MOSl/SW thin 


33710 




37070 


53636 


2005 Sep 23 


52.76 


S2 


EPIC-M0S2/SW thin 


33977 




37706 


53636 


2005 Sep 23 


52.77 


6369^ 


HRC-S/LETG 




7696 




53652 


2005 Oct 8 


26.26 


71779 


HRC-S/LETG 




2532 




53653 


2005 Oct 9 


8.04 


1086 S3 


EPIC-pn/FF thin 


167725 




117828 


53687 


2005 Nov 12 


37.84 
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Table 1. - Continued. 



Obs. Id. 


Instrument/setup 


Counts 




Counts 


MJD 


Start Date 


Effective 






(120-400 e^ 


V) (400-1000 cV) 






exposure 












[days] 




[ks] 


SI 


EPIC-MOSl/SW thin 


25265 




27091 


53687 


2005 Nov 12 


39.46 


S2 


EPIC-M0S2/SW thin 


25045 




27580 


53687 


2005 Nov 12 


39.47 


7243^° 


HRC-S/LETG 




5931 




53718 


2005 Dec 14 


17.18 


724410 


HRC-S/LETG 




4927 




53718 


2005 Dec 15 


16.29 


724510 


HRC-S/LETG 




5249 




53718 


2005 Dec 16 


17.19 


558410 


HRC-S/LETG 




5321 




53718 


2005 Dec 17 


14.19 


7251 


HRC-S/LETG 




^787 




53775 


2006 Sep 9 


10.65 


1181 SI 


EPIC-pn/FF thin 


89015 




59524 


53877 


2006 May 22 


20.04 


S2 


EPIC-MOSl/SW thin 


13652 




14219 


53877 


2006 May 22 


21.66 


S3 


EPIC-M0S2/SW thin 


15045 




15048 


53877 


2006 May 22 


21.66 


1265 SI 


EPIC-pn/FF thin 


89852 




60883 


54044 


2006 Nov 5 


20.04 


S2 


EPIC-MOSl/FF thin 


12414 




12296 


54044 


2006 Nov 5 


21.61 


S3 


EPIC-M0S2/FF thin 


15291 




14653 


54044 


2006 Nov 5 


21.62 


1356 SI 


EPIC-pn/FF thin 


89411 




57483 


54226 


2007 May 5 


20.04 


S2 


EPIC-MOSl/FF thin 


16690 




16685 


54226 


2007 May 5 


21.61 


S3 


EPIC-M0S2/FF thin 


17147 




16818 


54226 


2007 May 5 


21.62 


1454 SI 


EPIC-pn/FF thin 


102917 




64833 


54421 


2007 Nov 17 


23.06 


S2 


EPIC-MOSl/FF thin 


18674 




18197 


54421 


2007 Nov 17 


24.62 


S3 


EPIC-M0S2/FF thin 


19029 




18416 


54421 


2007 Nov 17 


24.62 


10861" 


HRC-S/LETG 




42AQ 




54851 


2009 Jan 20 


11.91 


10700" 


HRC-S/LETG 




14053 




54876 


2009 Feb 14 


21.82 


1700 S3 


EPIC-pn/FF thin 


37085 




30806 


54913 


2009 Mar 21 


10.84 


U2 


EPIC-MOSl/FF thin 


11481 




11876 


54912 


2009 Mar 21 


17.69 


U2 


EPIC-M0S2/FF thin 


11936 




11556 


54912 


2009 Mar 21 


17.69 


10701^2 


HRC-S/LETG 




15647 




55086 


2009 Sep 11 


33.17 


1792 S3 


EPIC-pn/FF thin 


65407 




48105 


55096 


2009 Sep 22 


17.90 


SI 


EPIC-MOSl/SW thin 


12524 




11919 


55096 


2009 Sep 22 


19.22 


S2 


EPIC-M0S2/SW thin 


12724 




11993 


55096 


2009 Sep 22 


19.24 



Note: Observation 7251, in italic, is not used in this work (see Sec. 3). The identifier marks the merged data sets: 
observations with the same tag are merged into one data set. Counts for ROSAT and Chandra HRC were not divided 
into soft (120-400 eV) and hard (400-1000 eV) band (see Sec. 3), hence we list the total number of counts. The count 
numbers for the HRC observations are comparable to those listed in KvK05, but are a factor 100 less than those listed 
in vK07 (probably due to missprint in vK07) although the data handling in all three works (KvK05, vK07 and this 
work) is comparable. The count numb ers of all other obse rvat ions are in good a greement. For details on ROSAT and 
XMM-Newton observations we refer to [Cropper et al.||2004 and Hohle et al. 2009 (and Sec. 2, 3 this work), respectively. 
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Hohle et al.: Timing solution of RX J0720.4-3125 



Table 2. Individual periods and TOAs of the observations of RX J0720.4-312f 


J as listed in Tab. [T] 




Obs. Id. 


Instrument/setup 


period 


period 


TOA 


TOA 






(120-400 eV) 


(400-1000 eV) 


(120-400 eV) 


(400-1000 eV) 






[s] 


[s] 


[days] 


[days' 


rp300338n00 


PSPC 


8.39120(44) 


49,257.2547153(16) 


49,257.2547298(21) 


rh300508n00 


HRI 


8.3852(63) 


50,198.6873383(25) 


50,198.6873509(35) 


rhlSOlOOnOO 


HRI 


8.3443(11) 


50,210.5562791(18) 


50,210.5562812(17) 


rh300508n01 


HRI 


8.4902(98) 


50,353.9975633(27) 


50,353.9975702(20) 


rh400884n00 


HRI 


8.391130(50) 


50,391.3004644(14) 


50,391.3004729(18) 


h400944n00 


HRI 


8.3921(10) 


50,925.6878393(11) 


50,925.6878472(21) 


368^ 


HRC-S/LETG 


8.39063(49) 


51,577.0395641(21) 


51,577.0395693(16) 


745I 


HRC-S/LETG 










3691 


HRC-S/LETG 










0078 S3 


EPIC-pn/FF thin 


8.391113(18) 


8.391085(35) 


51,677.44324000(30) 


51,677.44323979(69) 


SI 


EPIC-MOSl/FF thin 


8.391090(45) 


8.391090(55) 


51,677.47179274(83) 


51,677.47179286(78) 


S2 


EPIC-M0S2/SW thin 


8.391113(23) 


8.391033(53) 


51,677.47179234(54) 


51,677.4717923(13) 


0175 S3 


EPIC-pn/FF med 


8.391268(58) 


8.39124(12) 


51,869.95710586(60) 


51,869.95710803(84) 


S7 


EPIC-MOSl/LW med 


8.39122(22) 


8.39107(32) 


51,869.9949836(13) 


51,869.9949874(27) 


27742 


ACIS-CC 


8.391093(13) 


8.391133(13) 


52,248.6767290(15) 


52,248.67672069(69) 


27732 


ACIS-CC 










2771^ 


ACIS-CC 










27722 


ACIS-CC 










0533 S3 


EPIC-pn/FF thin 


8.391143(43) 


8.391098(73) 


52,584.92605120(69) 


52,584.9260530(10) 


SI 


EPIC-MOSl/FF thin 


8.39136(12) 


8.39088(16) 


52,584.91993465(11) 


52,584.9200294(15) 


S2 


EPIC-M0S2/FF thin 


8.39122(12) 


8.39104(20) 


52,584.91993465(95) 


52,584.9199325(16) 


0534 S3 


EPIC-pn/FF thin 


8.391225(40) 


8.391317(63) 


52,587.00129529(39) 


52,587.00129506(97) 


SI 


EPIC-MOSl/FF thin 


8.39119(11) 


8.39116(18) 


52,586.9952779(18) 


52,586.9952766(29) 


S2 


EPIC-M0S2/FF thin 


8.39137(13) 


8.39137(23) 


52,586.99527264(93) 


52,586.9952735(14) 


0622 U2 


EPIC-pn/SW thick 


8.391133(18) 


8.391138(23) 


52,761.99514720(82) 


52,761.99514812(64) 


S5 


EPIC-MOSl/FF med 


8.39130(16) 


8.39100(20) 


52,762.2413478(13) 


52,762.2413532(20) 


S3 


EPIC-MOSl/FF thin 


8.39101(13) 


8.39100(17) 


52,761.8668566(19) 


52,761.8668516(28) 


S6 


EPIC-M0S2/FF med 


8.39100(15) 


8.39126(15) 


52,762.2413439(17) 


52,762.24144068(14) 


S4 


EPIC-M0S2/FF thin 


8.390950(95) 


8.39118(21) 


52,761.8669520(11) 


52,761.8669491(20) 


0711 S8 


EPIC-pn/SW med 


8.391082(48) 


8.391188(82) 


52,940.11635982(97) 


52,940.11627407(69) 


S3 


EPIC-MOSl/FF thin 


8.39125(31) 


8.39077(52) 


52,939.99012307(84) 


52,939.9900211(20) 


S5 


EPIC-MOSl/LW thin 


8.39087(24) 


8.39059(41) 


52,940.1679353(12) 


52,940.1679433(29) 


S4 


EPIC-M0S2/FF thin 


8.39169(32) 


8.39051(49) 


52,939.9900070(17) 


52,939.9900222(24) 


S6 


EPIC-M0S2/LW thin 


8.39091(22) 


8.39147(26) 


52,940.1679318(13) 


52,940.1679407(14) 


46663 


ACIS-CC 


8.3911142(23) 


8.3911216(28) 


53,016.68131016(91) 


53,016.6813159(12) 


46673 


ACIS-CC 










4668^ 


ACIS-CC 










46693 


ACIS-CC 










5305^ 


HRC-S/LETG 


8.39109(25) 


53,062.4157092(15) 


53,062.4157091(15) 


0815 SI 


EPIC-pn/FF thin 


8.391105(35) 


8.391143(58) 


53,147.68119387(43) 


53,147.6812041(15) 


S2 


EPIC-MOSl/FF thin 


8.391108(68) 


8.39111(15) 


53,147.6883801(11) 


53,147.6883883(31) 


S3 


EPIC-M0S2/FF thin 


8.391150(60) 


8.39120(13) 


53,147.68838022(93) 


53,147.6883911(21) 


4670^ 


ACIS-CC 


8.3911146(15) 


8.3911165(16) 


53,230.5756209(11) 


53,230.5756243(19) 


4671^ 


ACIS-CC 










4672^ 


ACIS-CC 










4673^ 


ACIS-CC 










5581*^ 


HRC-S/LETG 


8.391010(95) 


53,393.6674769(19) 


53,393.6674775(23) 


0986 S3 


EPIC-pn/FF thin 


8.391138(23) 


8.391140(30) 


53,488.67561302(44) 


53,488.67561797(84) 


SI 


EPIC-MOSl/SW thin 


8.391095(45) 


8.391008(48) 


53,488.6695937(14) 


53,488.6695969(13) 


S2 


EPIC-M0S2/SW thin 


8.391183(38) 


8.391018(53) 


53,488.6696892(13) 


53,488.6696931(14) 


5582^ 


HRC-S/LETG 


8.391 


138(63) 


53,522.9398395(35) 


53,522.9398385(29) 


6364^ 


HRC-S/LETG 


8.39129(25) 


53,610.0881154(30) 


53,610.0881131(20) 


1060 S3 


EPIC-pn/FF thin 


8.391120(20) 


8.391082(33) 


53,636.30015746(23) 


53,636.3001616(10) 


SI 


EPIC-MOSl/SW thin 


8.391095(40) 


8.391130(50) 


53,636.29413571(70) 


53,636.2941416(13) 


S2 


EPIC-M0S2/SW thin 


8.391103(38) 


8.391163(58) 


53,636.29413754(96) 


53,636.29413962(83) 


6369^ 


HRC-S/LETG 


8.391 


10(11) 


53,652.2601916(17) 


53,652.2601885(15) 


71779 


HRC-S/LETG 










1086 S3 


EPIC-pn/FF thin 


8.391125(25) 


8.391072(38) 


53,687.17180535(53) 


53,687.17180817(92) 
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Table 2. - Continued. 



Obs. Id. 


Instrument /setup 


period 


period J 


TOA 


TOA 






(120-400 eV) 


(400-1000 eV) 


(120-400 eV) 


(400-1000 eV) 






¥\ 


[s] 


[days] 


[days] 


SI 


EPIC-MOSl/SW thin 


8.391003(63) 


8.391085(80) 


53,687.1657835(12) 


53,687.16588680(52) 


S2 


EPIC-M0S2/SW thin 


8.391200(55) 


8.390998(78) 


53,687.1658825(10) 


53,687.1658857(12) 


724310 


HRC-S/LETG 


8.391115(10) 


53,720.0243035(13) 


53,720.0243020(12) 


724410 


HRC-S/LETG 










7245 10 


HRC-S/LETG 










558410 


HRC-S/LETG 










7251 


HRC-S/LETG 


8.39212(88) 


53,775.3509131(17) 


53,775.35091144(12) 


1181 SI 


EPIC-pn/FF thin 


8.391090(70) 


8.39120(11) 


53,877.32857154(71) 


53,877.3285752(11) 


S2 


EPIC-MOSl/SW thin 


8.39137(17) 


8.39101(20) 


53,877.3225489(15) 


53,877.32255272(65) 


S3 


EPIC-M0S2/SW thin 


8.39131(13) 


8.39116(18) 


53,877.3225491(11) 


53,877.3225524(16) 


1265 SI 


EPIC-pn/FF thin 


8.391075(65) 


8.39104(11) 


54,044.60537610(44) 


54,044.6053791(11) 


S2 


EPIC-MOSl/FF thin 


8.39117(16) 


8.39074(17) 


54,044.5990605(12) 


54,044.5990664(13) 


S3 


EPIC-M0S2/FF thin 


8.39099(13) 


8.39096(19) 


54,044.59906220(62) 


54,044.5990665(18) 


1356 SI 


EPIC-pn/FF thin 


8.391150(65) 


8.391113(98) 


54,225.84098865(82) 


54,225.84099190(98) 


S2 


EPIC-MOSl/FF thin 


8.39107(21) 


8.39066(27) 


54,225.8345794(18) 


54,225.83457862(99) 


S3 


EPIC-M0S2/FF thin 


8.39103(23) 


8.39096(22) 


54,225.8346806(13) 


54,225.8346782(22) 


1454 SI 


EPIC-pn/FF thin 


8.391130(50) 


8.391085(85) 


54,421.36989356(53) 


54,421.36989479(60) 


S2 


EPIC-MOSl/FF thin 


8.39090(16) 


8.39109(17) 


54,421.3635797(17) 


54,421.3635840(15) 


S3 


EPIC-M0S2/FF thin 


8.39131(22) 


8.39127(22) 


54,421.3635875(20) 


54,421.3635849(18) 


10861" 


HRC-S/LETG 


8.39129(20) 


54,863.7916775(17) 


54,863.7916751(12) 


10700" 


HRC-S/LETG 










1700 S3 


EPIC-pn/FF thin 


8.39108(99) 


8.3912(12) 


54,911.7856154(16) 


54,911.78046975(81) 


U2 


EPIC-MOSl/FF thin 


8.39097(26) 


8.39136(29) 


54,911.7440458(30) 


54,911.7440479(17) 


U2 


EPIC-M0S2/FF thin 


8.39103(25) 


8.39118(35) 


54,911.7440524(13) 


54,911.7440483(15) 


10701^2 


HRC-S/LETG 




8.39119(20) 


55,085.6684206(14) 


55,085.6684187(15) 


1792 S3 


EPIC-pn/FF thin 


8.39102(12) 


8.39101(15) 


55,096.30446151(69) 


55,096.3044632(11) 


SI 


EPIC-MOSl/SW thin 


8.39097(20) 


8.39120(20) 


55,096.2966932(13) 


55,096.2966943(13) 


S2 


EPIC-M0S2/SW thin 


8.39104(15) 


8.39102(24) 


55,096.2967896(14) 


55,096.2967885(11) 



Note: The individual periods for ROSAT and Chandra HRC were not divided into soft (120-400 eV) and hard (400- 
1000 eV) band (see Sec. 3), hence we hst the individual periods derived from all counts. We hst the TOAs (definition 
equal to that in KvK05 and vK07) of the soft band (for the "all data" solution of the soft ban d, see Tab. [3| and o f the 
hard band (for the "all data" solution of the hard band, see Tab. [3]). The errors of the periods ( Ransom et "al.|2002 1 and 
the TOAs (derived from the fitted light curves, see text) are given in parenthesis and denote la confidence level. 
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